Primordial non-Gaussianity in the Large Scale Structure of the Universe 



O 
(N 

(N 



Vincent Desjacques^' * and Uros Seljak^' ^' 

'institute for Theoretical Physics, University of Zurich, 8057 Zurich, Switzerland 
''Physics Department, Astronomy Department and Lawrence Berkeley National Laboratory, 
University of California, Berkeley, California 94720, USA 
(Dated: June 25, 2010) 

Primordial non-Gaussianity is a potentially powerful discriminant of the physical mechanisms that 
generated the cosmological fluctuations observed today. Any detection of significant non-Gaussianity 
would thus have profound implications for our understanding of cosmic structure formation. The 
large scale mass distribution in the Universe is a sensitive probe of the nature of initial conditions. 
Recent theoretical progress together with rapid developments in observational techniques will enable 
us to critically confront predictions of infiationary scenarios and set constraints as competitive as 
those from the Cosmic Microwave Background. In this paper, we review past and current efforts in 
the search for primordial non-Gaussianity in the large scale structure of the Universe. 
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I. INTRODUCTION 

In generic inflationary models based on the slow roll 
of a scalar field, primordial curvature perturbations are 
produced by the inflaton field as it slowly rolls down its 
potential [1-4]. Most of these scenarios predict a nearly 
scale-invariant spectrum of adiabatic curvature fluctu- 
ations, a relatively small amount of gravity waves and 
tiny deviations from Gaussianity in the primeval distri- 
bution of curvature perturbations [5-7]. Although the 
latest measurements of the cosmic microwave background 
(CMB) anisotropics favor a slightly red power spectrum 
[8], no significant detection of a B-mode or some level 
of primordial non-Gaussianity (NG) has thus far been 
reported from CMB observations. 

While the presence of a B-mode can only be tested 
with CMB measurements, primordial deviations from 
Gaussianity can leave a detectable signature in the dis- 
tribution of CMB anisotropics and in the large scale 
structure (LSS) of the Universe. Until recently, it was 
widely accepted that measurement of the CMB furnish 
the best probe of primordial non-Gaussianity [9] (see, 
e.g., the recent review by E. Komatsu on primordial 
non-Gaussianity in the CMB [10]). However, these con- 
clusions did not take into account the anomalous scale- 
dependence of the galaxy power spectrum and bispec- 
trum arising from primordial NG of the local type 
[11, 12]. These theoretical results, together with rapid de- 
velopments in observational techniques, will provide large 
amount of LSS data to critically confront predictions of 
non-Gaussian models. In particular, galaxy clustering 
should provide independent constraints on the magnitude 
of primordial non-Gaussianity as competitive as those 
from the CMB and, in the long run, may even give the 
best constraints. 



The purpose of this work is to provide an overview of 
the search for a primordial non-Gaussian signal in the 
large scale structure. We will begin by briefly summa- 
rizing how non-Gaussianity arises in inflationary models 
(§11). Next, we will discuss the impact of primordial non- 
Gaussianity on the mass distribution in the low redshift 
Universe (§111). The main body of this review is §IV, 
where we describe in detail an number of methods ex- 
ploiting the abundance and clustering properties of ob- 
served tracers of the LSS to constrain the amount of ini- 
tial non-Gaussianity. We conclude with a discussion of 
present and forecasted constraints achievable with LSS 
surveys (§V). 



II. MODELS AND OBSERVABLES 

Because they assume i) a single dynamical fleld (the 
inflaton) ii) canonical kinetic energy terms (i.e. per- 
turbations propagate at the speed of light) iii) slow roll 
(i.e. the timescale over which the inflaton field changes 
is much larger than the Hubble rate) iv) an initial vac- 
uum state, single-field slow-roll models lead to a small 
level of primordial non-Gaussianity [6, 7, 13]. The lowest 
order statistics sensitive to non-Gaussian features in the 
initial distribution of scalar perturbations 3>(x) (we shall 
adopt the standard CMB convention in which <I>(x) is the 
Bardeen's curvature perturbation in the matter era) is 
the 3-point function or bispectrum i3$(ki, k2, ks), which 
is a function of any triangle ki -f k2 -t- ks = (as follows 
from statistical homogeneity which we assume through- 
out this paper). It has been recently shown that, in the 
squeezed limit <^ ki m k2-, the bispectrum of any 
single-field slow-roll inflationary model asymptotes to the 
local shape (defined in Eq. 3) [14-16]. The corresponding 
nonlinear parameter predicted by these models is 
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where Ug is the tilt or spectral index of the power spec- 
trum P$(fc), which is accurately measured to be rig « 
0.960 ± 0.013 [8]. Therefore, any robust measurement of 
/nl well above this level would thus rule out single-field 
slow-roll inflation. 



A. The shape of the primordial bispectrum 

Large, potentially detectable amount of Gaussianity 
can be produced when at least one of the assumptions 
listed above is violated, i.e. by multiple scalar flelds 
[17, 18], nonlinearities in the relation between the pri- 
mordial scalar pcrtiubations and the inflaton field [7, 13], 
interactions of scalar fields [19], a modified dispersion re- 
lation or a departure from the adiabatic Bunch-Davies 
ground state [20], Generation of a large non-Gaussian 
signal is also expected at reheating [21] and in the ekpy- 
rotic scenario [22]. Each of these physical mechanisms 
leaves a distinct signature in the primordial 3-point func- 
tion _B$(ki, k2, ka), a measurement of which would thus 
provide a wealth of information about the physics of pri- 
mordial fluctuations. Although the conflguration shape 
of the primordial bispectrum can be extremely complex 
in some models, there are broadly three classes of shape 
characterizing the local, equilateral and folded type of 
primordial non-Gaussianity [23, 24]. The magnitude of 
each template "X" is controlled by a dimensionless non- 
linear parameter /^^ which we seek to constrain us- 
ing CMB or LSS observations (instead of attempting a 
model- independent measurement of 

Any non-Gaussianity generated outside the horizon in- 
duces a 3-point function that is peaked on squeezed or 
collapsed triangles for realistic values of the scalar spec- 
tral index. The resulting non-Gaussianity depends only 
on the local value of the Bardeen's curvature potential, 
and can thus be conveniently parameterized up to third 
order by [7, 9, 13, 25] 



$(x) 



(x)+/^-02(x) 



(2) 



where 0(x) is an isotropic Gaussian random fleld and 
/nL' 3nl ^'^e dimensionless, phenomenological parame- 
ters. Since curvature perturbations are of magnitude 
O(10~^), the cubic order correction should always be 
negligibly small compared to the quadratic one when 
^(./nl) ^ ^(.9nl)- However, this condition is not sat- 
isfied by some multifield inflationary models such as the 
curvaton scenario, in which a large (jr|^£ and a small 
can be simultaneously produced [18] . The quadratic term 
generates the 3-point function at leading order, 

Sl,-(ki,k2,k3) = 2f^°l [P^{h)P^{k2) + (eye.)] , (3) 

where (eye.) denotes all cyclic permutations of the in- 
dices and P^{k) is the power spectrum of the Gaussian 
part <?i>(x) of the Bardeen potential. The cubic-order 
terms generates a trispectrum T$(ki,k2,k3,k4) at lead- 
ing order. 



Equilateral type of non-Gaussianity, which arises in 
inflationary models with higher-derivative operators such 
as the DBI model, is well described by the factorizable 
form [26] 



B^''(ki,k2,k3) = 6/^1 -(P^(fci)P^(fc2) + (eye.)) 
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(4) 



+ {P;'\k,)P^''\k2)P4k3) + (perm.)) 

It can be easily checked that the signal is largest in the 
equilateral configurations ki « ^2 « fcs, and suppressed 
in the squeezed limit k^ <C fci ~ ^2. Note that, in single- 
field slow-roll infiation, the 3-point function is a linear 
combination of the local and equilateral shape [14] . 

As a third template, we consider the folded or flattened 
shape which is maximized for A;2 ~ fca ~ A;i/2 [27] 



B4°Hki,k2,k3) = 6f^t (P4fci)P4fc2) + (eye.)) 



+3{P^{ki)P^{k2)PM) 
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(5) 



- {P'^"{k^)Pl"{k2)PM + (perm.)) 

and approximate the non-Gaussianity due to modifica- 
tion of the initial Bunch-Davies vacuum in canonical sin- 
gle field action (the actual 3-point function is not factor- 
izable). As in the previous example, is suppressed 
in the squeezed configurations. Unlike B'^ however, 
the folded shape induces a scale-dependent bias at large 
scales (see §IVC). 



B. Statistics of the linear mass density fleld 

The Bardeen's curvature potential $(x) is related to 
the linear density perturbation (5o(k, 0) at redshift z 
through 



5o{\<i,z)=M{k,z)^{\i) 



where 



M{k,z) = 



2 k^T(k)D{z) 

3 Om^fo 



(6) 



(7) 



Here, T{k) is the matter transfer function normalized to 
unity as fc — ^ 0, fim is the present-day matter density 
and D{z) is the linear growth rate normalized to 1 -|- 
n-point correlators of the linear mass density field can 
thus be written as 

((5o(k) • • • 5o(k„)) = (jlMih)^ ($(ki) • • • $(k„)) . 

(8) 

Smoothing inevitably arises when comparing observa- 
tions of the large scale structure with theoretical pre- 
dictions from, e.g., perturbation theory (PT) which are 
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valid only in the weakly nonlinear regime [28] , or from the 
spherical collapse model which ignores the strongly non- 
linear internal dynamics of the collapsing regions [29, 30]. 
For this reason, we introduce the smoothed linear density 
field 

SniK z) = M{k, 2)l^fl(A:)$(k) = Mnik, z)$(k) , (9) 

where Wii{k) is a (spherically symmetric) window func- 
tion of characteristic radius R or mass scale M that 
smoothes out the small-scale nonlinear fluctuations. We 
will assume a top-hat filter throughout. Furthermore, 
since M and R are equivalent variables, we shall indis- 
tinctly use the notation Sr and Sm in what follows. 



C. Topological defects models 

In addition to inflationary scenarios, there is a whole 
class of models, known as topological defcict models, in 
which cosmological fluctuations are sourced by an inho- 
mogeneously distributed component which contributes 
a small fraction of the total energy momentum tensor 
[31, 32]. The density field is obtained as the convolu- 
tion of a discrete set of points with a specific density 
proflle. Defects are deeply rooted in particle physics as 
they are expected to form at a phase transition. Since 
the early Universe may have plausibly undergone several 
phase transitions, it is rather unlikely that no defects at 
all were formed. Furthermore, high redshift tracers of the 
LSS may be superior to CMB at finding non-Gaussianity 
sourced by topological defects [33]. However, CMB ob- 
servations already provide stringent limits on the energy 
density of a defect component [8] , so we shall only mini- 
mally discuss the imprint of these scenarios in the large 
scale structure. Phenomenological defect models are for 
instance 

5(x) = </.(x) + aNL (<A'(x) - (</>')) (10) 

in which the initial matter density 5(x) (rather than the 
curvature perturbation $(x)) contains a term propor- 
tional to the square of a Gaussian scalar field (f){x) [9], or 
the model (also known as isocurvature CDM model) 
in which ^(x) oc (?i^(x) [34]. 



III. EVOLUTION OF THE MATTER DENSITY 
FIELD WITH PRIMORDIAL NG 

In this Section, we summarize a number of results rel- 
ative to the effect of primordial NG on the mass density 
field. These will be useful to understand the complexity 
that arises when considering biased tracers of the density 
field (see §IV). 



A. Setting up non-Gaussian initial conditions 

Investigating the impact of non-Gaussian initial condi- 
tions (ICs) on the large scale structure traced by galaxies 
etc. requires simulations large enough so that many long 
wavelength modes are sampled. At the same time, the 
simulations should resolve dark matter halos hosting lu- 
minous red galaxies (LRGs) or quasars (QSOs), so that 
one can construct halo samples whose statistical proper- 
ties mimic as closely as possible those of the real data. 
This favors the utilization of pure N-body simulations for 
which a large dynamical range can be achieved. 

The evolution of the matter density field with primor- 
dial non-Gaussianity has been studied in series of large 
cosmological N-body simulations seeded with Gaussian 
and non-Gaussian initial conditions, see e.g. [12, 35-44]. 
For generic non-Gaussian (scalar) random fields, we face 
the problem of setting up numerical simulations with a 
prescribed correlation structure [45]. While an imple- 
mentation of the equilateral and folded bispectrum shape 
requires the calculation of several computationally de- 
manding convolutions, the operation is straightforward 
for primordial NG described by a local mapping such as 
the or the model. In the latter case, the local 
transformation $ = + /nl?^^ performed before mul- 
tiplication by the matter transfer fimction T{k) (com- 
puted with publicly available Boltzmann codes [46, 47]). 
The (dimensionless) power spectrum of the Gaussian 
part (f>{x) of the Bardeen potential is the usual power- 
law A2(fc) = fc3p^(fc)/(2^2) ^ A4k/kor^-\ Un- 
less otherwise stated, we shall assume a normalization 
= 7.96 X 10-1° at the pivot point fco = 0.02Mpc-^ To 
date, essentially all numerical studies of structure forma- 
tion with inflationary non-Gaussianity have implemented 
the local shape solely, so we will focus on this model in 
what follows. 

Non-Gaussian corrections to the primordial curvature 
perturbations can renormalize the input (unrenormal- 
ized) power spectrum of fluctuations used to seed the sim- 
ulations [48]. For the local fj^l model with \f^l\ < 100, 
renormalization effects are unlikely to be noticeable due 
to the limited dynamical range of current cosmological 
simulations. However, they can be significant, for ex- 
ample, in simulations of a local cubic coupling g^l^<j)^ 
with a large primordial trispectrum [49]. The cubic or- 
der term g^f^'^ renormalizes the amplitiidc of the 
power spectrum of initial curvature perturbations to 
A^^ + 65fjj£((^2), where 




For scale invariant initial conditions, ((/)^) has a logarith- 
mic divergence at large and small scales. In practice how- 
ever, the finite box size and the resolution of the simu- 
lations naturally furnish a low- and high-A; cutoff. The 
effective correction to the amplitude of density fluctua- 
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FIG. 1: Skewness cr Sg^^T?) of the smoothed density field in 
unit of /^L for the local, equilateral and folded bispectrum 
shape. The skewness for the equilateral and folded templates 
is a factor of ~ 3 smaller than in the local model. In any 
case, this implies that \(JrSz(R)\ ^ 1 on the scales probed 
by the large scale structure for realistic values of the non- 
linear coupling parameter, |/nlI ^ 100. The shaded regions 
approximately indicate the range of scales probed by various 
LSS tracers. For the galaxy power spectrum and bispectrum, 
the upper limit sensitively depends upon the surveyed volume. 
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tions (5(78 in the mo(iel thus is 

where N is the number of mesh points along one dimen- 
sion and i is the simulation box length. For ^I^l ~ 
L = 1.5 /i^^Gpc and N = 1024 for instance, we obtain 
Sas « 0.015. 

To generate the initial particle distribution, the Zel- 
dovich approximation is commonly used instead of the 
exact gravitational dynamics. This effectively corre- 
sponds to starting from non-Gaussian initial conditions 
[50]. Since the transition to the true gravitational dy- 
namics proceeds rather gradually [51], one should ensure 
that the initial expansion factor is much smaller than 
that of the outputs analyzed. Alternatively, it is possi- 
ble to generate more accurate ICs based on second-order 
Lagrangian perturbation theory (2LPT) [52]. At fixed 
initial expansion factor, they reduce transients such that 
the true dynamics is recovered more rapidly [53] . 



B. Mass density probability distribution 

In the absence of primordial NG, the probability dis- 
tribution function (PDF) of the initial smoothed density 



field (the probability that a randomly placed cell of vol- 
ume V has some specific density) is Gaussian. Namely, 
all normalized or reduced smoothed cumulants Sj of or- 
der J > 3 are zero. An obvious signature of primordial 
NG would thus be an initially non-vanishing skewness 
^3 = {S'j,)e/{5ir or kurtosis ^4 = (4)^/(^1^)' " V{Sl) 
[37, 54, 55]. Here, the subscript c denotes the connected 
piece of the n-point moment that cannot be simplified 
into a sum over products of lower order moments. At 
third order for instance, the cumulant of the smoothed 
density field is an integral of the 3-point function. 



where 



d^ki fd^k2 fd^k. 



(27r)3 J (27r)3 J (27r) 



Bfl(ki,k2,k3,z) , 

(13) 



Bij(ki,k2,k3,z) = MR{ki,z)MR{k2,z)MR{k3,z) 

X B$(ki,k2,k3) (14) 

is the bispectrum of the smoothed linear density fluc- 
tuations at redshift z. Note that, owing to S3{R, z) oc 
D{z)~^, the product aS3{R) does not depend on redshift. 
Over the range of scale accessible to LSS observations, 
aS3{R) is a monotonic decreasing function of R that is 
of magnitude ^ 10^'* for the local, equilateral and folded 
templates discussed above (Figure 1). Strictly speaking, 
all reduced moments should be specified to fully charac- 
terize the density PDF, but a reasonable description of 
the density distribution can be achieved with moments 
up to the fourth order. 

Numerical and analytic studies generally find that a 
density PDF initially skewed towards positive values pro- 
duces more overdense regions while a negatively skewed 
distribution produces larger voids. Gravitational insta- 
bilities also generate a positive skewness in the density 
field, refiecting the fact that the evolved density distribu- 
tion exhibits an extended tail towards large overdensities 
[50, 56-60]. This gravitationally-induced signal eventu- 
ally dominates the primordial contribution such that, at 
fixed normalization amplitude, non-Gaussian models de- 
viate more strongly from the Gaussian paradigm at high 
redshift. The time evolution of the normalized cumu- 
lants Sj can be worked out for generic Gaussian and 
non-Gaussian ICs using, e.g., PT or the spherical collapse 
approximation. For Gaussian ICs, PT predicts the nor- 
malized cumulants to be time-independent to lowest non- 
vanishing order, with a skewness S3 ~ 34/7, whereas for 
non-Gaussian ICs, the linear contribution to the cumu- 
lants decays as Sj{R,z) = Sj{R,oo)/D-^-^{z) [61, 62]. 

The persistence of the primordial hierarchical ampli- 
tude 5*7 (i?, 00) sensitively depends upon the magnitude 
of Sn, N > J, relative to unity. For example, an ini- 
tially large non- vanishing kurtosis could source skewness 
with a time-dependence and amplitude similar to that 
induced by nonlinear gravitational evolution [61]. Al- 
though there is an infinite class of non-Gaussian models, 
we can broadly divide them into weakly and strongly non- 
Gaussian. In weak NG models, the primeval signal in the 
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normalized cumulants is rapidly obliterated by gravity- 
induced non-Gaussianity. This is the case of hierarchical 
scaling models where n-point correlation functions sat- 
isfy ^„ oc ^2 ~^ with ^2 ^ 1 at large scales. By contrast, 
strongly NG initial conditions dominate the evolution of 
the normalized cumulants. This occurs when the hierar- 
chy of correlation functions obeys the dimensional scaling 

^ ^2^^! which arises in the particular case of ini- 
tial conditions [63] or in defect models such as texture 
[38, 64, 65]. These scaling laws have been successfully 
confronted with numerical investigations of the evolution 
of cumulants [38, 39]. We note that the scaling of the 
contribution induced by gravity is, however, different for 
the kurtosis [66], suggesting that the latter is a better 
probe of the nature of initial conditions. 

Although gravitational clustering tends to erase the 
memory of initial conditions, numerical simulations of 
non-Gaussian initial conditions show that the occurrence 
of highly underdense and overdense regions is very sensi- 
tive to the presence of primordial NG. In fact, the imprint 
of primordial NG is best preserved in the negative tail 
of the PDF P{pr) of the evolved (and smoothed) den- 
sity field PR [41]. A satisfactory description of this effect 
can be obtained from an Edgeworth expansion of the ini- 
tial smoothed overdensity field [67]. At high densities 
PR ^ 1, the non-Gaussian modification approximately 

3/5 

scales as whereas, at low densities p;? ~ 0, the devi- 
ation is steeper, P^r^ ■ Taking into account the weak scale 
dependence of aSj,{R) further enhances this asymmetry. 
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FIG. 2: Non-Gaussian fractional correction /3s{k,z) — 
APg"^ {k, z) /Pp{k, z) to the matter power spectrum that orig- 
inates from primordial non-Gaussianity of the local type. Re- 
sults are shown at redshift z = 0, 0.5, 1 and 2 for f^l = -|-100 
(filled symbols) and /nl = —100 (empty symbols). The solid 
curves indicate the prediction from a 1-loop perturbative ex- 
pansion. 



C. Power spectrum and bispectrum 

Primordial non-Gaussianity also imprints a signature 
in Fourier space statistics of the matter density field. Pos- 
itive values of tend to increase the small scale mat- 
ter power spectrum Ps{k) [11, 41, 68] and the large scale 
matter bispectrum ^^(ki, k2, ka) [11, 69]. In the weakly 
nonlinear regime where 1-loop PT applies, the Fourier 
mode of the density field for growing-mode initial condi- 
tions reads [57, 70] 

J(k,z) = ^o(k,z) + j£qid\2SD{k- qi - ^2) 

X F2(qi,q2)5o(qi,2)(5o(q2,2:) . (15) 

The kernel f2(ki,k2) = 5/7+p{ki/k2+k2/ki)/2+2p'^/7, 
where p, is the cosine of the angle between ki and k2, 
describes the nonlinear 2nd order evolution of the density 
field. It is nearly independent of J7m and JIa and vanishes 
in the (squeezed) limit ki = — k2 as a consequence of 
the causality of gravitational instability. At 1-loop PT, 
Eq.(15) generates the mass power spectrum 

Ps{k, Z) = Pf{k, Z) + ApNG(fc^ ^) = ^) 

+ [P(22){k.z) + P^^s){k,z)\ +APf''ik,z) . (16) 



Here, Po(fc, z) is the linear matter power spectrum at 
redshift z, P(22) and ^"(13) are the standard one-loop con- 
tributions in the case of Gaussian ICs [70, 71], and 

APi^G z) = 2 1 ^F2 (q, k - q)Po (-k, q, k - q, z) 

(17) 

is the correction due to primordial NG [68]. This terms 
scales as cx D^{z), so the effect of non-Gaussianity is 
largest at low redshift. Moreover, because F2(ki,k2) 
vanishes in the squeezed limit, Eq.(17) is strongly sup- 
pressed at small wavenumbers, even in the local f^l 
model for which i3o(— k, q, k — q, z) is maximized in the 
same limit (i.e. |k| 0). For fj^£ - ©(lO^), the magni- 
tude of this correction is at a per cent level in the weakly 
nonlinear regime k < 0.1 hMpc~^ [42, 44, 72], in good 
agreement with the measurements (see Figure 2). Ex- 
tensions of the renormalization group description of dark 
matter clustering [73] to non-Gaussian initial density and 
velocity perturbations can improve the agreement up to 
wavenumbers k < 0.25 /iMpc~^ [74, 75]. 

It is also instructive to compare measurements of the 
matter bispectrum Bs{ki,k2,k3) with perturbative pre- 
dictions. To second order in PT, the matter bispectrum 
is the sum of a primordial contribution and two terms in- 
duced by gravitational instability [57, 76] (we will hence- 



6 



forth omit the cxphcit z-dependence for brevity), 

B5(ki,k2,k3) =So(ki,k2,k3) (18) 

+ [2F2(ki,k2)Po(fci)Po(A;2) + (eye.) 
F2(q,k3 - q) 



(27r)3 

X To(q, ka - q,ki,k2) + (eye.) 



where To(ki,k2,k3,k4) is the primordial trispeetrum of 
the density field. A similar expression can also be derived 
for the matter trispeetrum, which turns out to be less sen- 
sitive to gravitationally induced nonlinearities [77]. The 
reduced bispectrum Qs is conveniently defined as 



(53(ki,k2,k3) = 



B5(ki,k2,k3) 



Ps{h)Ps{k2) + cyclic 



(19) 



For Gaussian initial conditions, Q3 is independent of 
time and, at tree-level PT, is constant and equal to 
Q3{k,k,k) = 4/7 for equilateral configurations [57]. For 
general triangles moreover, it approximately retains this 
simple behavior, with a dependence on triangle shape 
through i^2(ki,k2) [11]. Figure 3 illustrates the effect 
of primordial NG of the local type on shape depen- 
dence of Q3 for a particular set of triangle configura- 
tions. As can be seen, the inclusion of 1-loop corrections 
greatly improve the agreement with the numerical data 
[78]. An important feature that is not apparent in Fig. 3 
is the fact that the primordial part to the reduced matter 
bispectrum scales as Q3 cx 1/Ai]i{k) for approximately 
equilateral triangles (and imder the assumption that /^l 
is scale-independent) [11]. This anomalous scaling con- 
siderably raises the ability of the matter bispectrum to 
constrain primordial NG of the local type. Unfortu- 
nately, neither the matter bispectrum nor the power spec- 
trum are directly observable with the large scale struc- 
ture of the Universe. Temperature anisotropics in the 
redshifted 21cm background from the pre-reionization 
epoch coiild in principle furnish a direct measurement 
of these quantities [79 81]. but foreground contamina- 
tion may severely hamper any detection. Weak lensing 
is another direct probe of the dark matter, although we 
can only observe it in projection along the line of sight 
[82]. 

As we will see shortly however, this large-scale anoma- 
lous scaling is also present in the bispectrum and power 
spectrum of observable tracers of the large scale struc- 
ture such as galaxies. It is this unique signature that will 
make future all-sky LSS surveys competitive with CMB 
experiments. 



D. Velocities 

Primordial non-Gaussianity also leaves a signature in 
the large-scale coherent bulk motions which, in the lin- 
ear regime, are directly related to the linear density field 



[56]. The various non-Gaussian models considered by 
[83] tend to have lower velocity dispersion and bulk flow 
than fiducial Gaussian model, regardless of the sign of 
the skewness. However, while the probability distribu- 
tion of velocity components is sensitive to primordial NG 
of the local type, in defect models it can be very close to 
Gaussian, even when the density field is strongly non- 
Gaussian, as a consequence of the central limit theorem 
[84, 85]. In this regards, correlation of velocity differ- 
ences could provide a better test of non-Gaussian initial 
conditions [86]. 

To measure peculiar velocities, one must subtract the 
Hubble flow from the observed redshift. This requires an 
estimate of the distance which is only available for nearby 
galaxies and clusters (although, e.g., the kinetic Sunyaev- 
Zel'dovich (kSZ) effect could be used to measure the bulk 
motions of distant galaxy clusters [87]). So far, measure- 
ments of the local galaxy density and velocity fields [88] 
as well as reconstruction of the initial density PDF from 
galaxy density and velocity data [89] are consistent with 
Gaussian initial conditions. 



IV. LSS PROBE OF PRIMORDIAL 
NON-GAUSSIANITY 

Discrete and continuous tracers of the large scale struc- 
ture such as galaxies, the Lya forest, the 21cm hydrogen 
line etc., provide a distorted image of the matter density 
field. In CDM cosmologies, galaxies form inside over- 
dense regions [90] and this introduce a bias between the 
mass and the galaxy distribution [91]. As a result, dis- 
tinct samples of galaxies trace the matter distribution 
differently, the most luminous galaxies preferentially re- 
siding in the most massive DM halos. This biasing effect, 
which concerns most tracers of the large scale structure, 
remains to be fully understood. Models of galaxy cluster- 
ing assume for instance that the galaxy biasing relation 
only depends on the local mass density, hut the actual 
mapping could be more complex [92, 93]. Because of bi- 
asing, tracers of the large scale structure will be aff'ected 
by primordial non-Gaussianity in a different way than the 
mass density field. In this Section, we describe a num- 
ber of methods exploiting the abundance and clustering 
properties of biased tracers to constrain the level of pri- 
mordial NG. We focus on galaxy clustering as it provides 
the tightest limits on primordial NG (see §V). 



A. Halo finding algorithm 

Locating groups of bound particles, or DM halos, in 
simulations is central to the methods described below. 
In practice, we aim at extracting halo catalogs with sta- 
tistical properties similar to those of observed galaxies or 
quasars. This, however, proves to be quite difficult be- 
cause the relation between observed galaxies an DM halos 
is somewhat uncertain. Furthermore, there is freedom at 
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FIG. 3: Reduced matter bispectrum Qa as a function of the 
angle Q between ki and k2 for a fixed k\ = 0.094 /iMpc~^ and 
hi = 1.5fci. The panels show ratios between the non-Gaussian 
and Gaussian Q3 for /n£ = +100 (top) and —100 (bottom). 
Dashed lines correspond to tree-level PT while continuous line 
indicate the 1-loop PT prediction. 



defining a halo mass. 

A important ingredient is the choice of the halo iden- 
tification algorithm. Halo finders can be broadly divided 
into two categories: Friends-of-Friends (FOF) finders [94] 
and spherical overdensity (SO) finders [95]. While the 
mass of a SO halo is defined by the radius at which the 
inner overdensity exceeds Avir(z) (typically ^ a few hun- 
dred times the background density ^(z)), the mass of a 
FOF halo is given by the number of particles within a 
linking length h from each other (6 ~ 0.15 — 0.2 in unit 
of mean interparticle distance) from each other. These 
definitions are somewhat arbitrary and may suit specific 
purposes only. In what follows, we shall mainly present 
results for SO halos as their mass estimate is more closely 
connected to the predictions of the spherical collapse 
model, on which most of the analytic formulae presented 
in this Section are based. The question of how the spher- 
ical overdensity masses can be mapped onto friends-of- 
friends masses remains a matter of debate (e.g. [96]). 
Clearly however, since the peak height i^(Af , z) depends 
on the halo mass M through the variance ctm (see below) , 
any systematic difference will be reflected in the value of 
V associated to a speciflc halo sample. As we will see 
shortly, this affects the size of the fractional deviation 
from the Gaussian mass function. 

Catalogs of mock galaxies with luminosities compara- 
ble to those of the targeted survey provide an additional 



layer of complication that can be used, among others, 
to assess the impact of observational errors on the non- 
Gaussian signal. However, most numerical studies of cos- 
mic structure formation with primordial NG have not yet 
addressed this level of sophistication, so we will discuss 
results based on statistics of dark matter halos only. 

B. Abundances of voids and bound objects 

It has long been recognized that departure from Gaus- 
sianity can significantly affect the abundance of highly 
biased tracers of the mass density field, as their frequency 
sensitively depends upon the tails of the initial density 
PDF [97-99]. The (extended) Press-Schcchter approach 
has been extensively applied to ascertain the effect of pri- 
mordial NG on the high mass tail of the mass function. 

1. Press-Schechter approach 

The Press-Schechter theory [100] and its extentions 
based on excursion sets [101-103] predict that the num- 
ber density n(M, z) of halos of mass M at redshift z is 
entirely specified by a multiplicity function f{v), 

where the peak height or significance i^(M, z) = 
dc{z)/aM is the typical amplitude of fiuctuations that 
produce those halos. Here and henceforth, ctm denotes 
the variance of the initial density field Sm smoothed on 
mass scale M (x and linearly extrapolated to present 
epoch, whereas Sc{z) « 1.68D{0)/ D{z) is the critical 
linear overdensity for (spherical) collapse at redshift z. 
In the standard Press-Schechter approach, n{M, z) is re- 
lated to the level excursion probability P{> dc,M) that 
the linear density contrast of a region of mass M exceeds 

<^c(^), 

where the last equality assumes Gaussian initial condi- 
tions. The factor of 2 is introduced to account for the 
contribution of low density regions embedded in over- 
densities at scale > M. In the extended Press-Schechter 
theory, Sm evolves with M and vf{v) is the probabil- 
ity that a trajectory is absorbed by the constant barrier 
5 = 5c (as is appropriate in the spherical collapse approx- 
imation) on mass scale M . In general, the exact form of 
f{v) depends on the barrier shape [104] and the filter 
shape [105]. Note also that jdv f{v) — 1, which ensures 
that all the mass is contained in halos. 

Despite the fact that the Press-Schechter mass func- 
tion overpredicts (underpredicts) the abundance of low 
(high) mass objects, it can be used to estimate the frac- 
tional deviation from Gaussianity. In this formalism. 
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FIG. 4: Fractional deviation from the Gaussian mass func- 
tion as a function of the peak height v = S^/a. Different 
symbols refer to different redshifts as indicated. The various 
curves are theoretical prediction at z = (see text). Ha- 
los were identified using a spherical overdensity (SO) finder 
with a redshift-dependent overdensity threshold A^iriz) (with 
Avir(z) increasing from ~ 200 at high redshift to attain ~ 350 
at 2: = 0). Error bars denote Poisson errors. For illustra- 
tion, M = 10^^ Mq/H corresponds to a significance v = 3.2, 
5.2, 7.7 at redshift z = 0, 1 and 2, respectively. Similarly, 
M = 10^* Mq/H and lO" Mq/H correspond to u = 1.9, 3, 
4.5 and 1.2, 1.9, 2.9 respectively. 



the non-Gaussian fractional correction to the multiplicity 
function is i?(t^,/^L) = fi'^JD/fi'^^O) = {dP/dM){> 
6c,MJ^^)/{dP/dM){> (5c,Af,0), which is readily com- 
puted once the non-Gaussian density PDF P{6m) is 
known. In the simple extensions proposed by [106] and 
[107], P{5m) is expressed as the inverse transform of a 
cumulant generating function. In [107], the saddle-point 
technique is applied directly to P{6m)- The resulting 
Edgeworth expansion is then used to obtain P(> Sc, M). 
Neglecting cumulants beyond the skewness, one obtain 
(we momentarily drop the subscript M for convenience) 



1 



ld{aS^) 



& d\a.v \ V 



(22) 



after integration over regions above bj^z). In [106], it is 
the level excursion probability P(> (Sc, M) that is calcu- 
lated within the saddle-point approximation. This ap- 
proximation better asymptotes to the exact large mass 
tail, which exponentially deviates from the Gaussian tail. 
To enforce the normalization of the resulting mass func- 
tion, one may define i^* — St, /a with 5* = 5c\/l — S^Sc/'i, 
and use [106, 108] 



2 n^^{M,z) dlnM 



din Vi, 



(23) 



The fractional deviation from the Gaussian mass function 
then becomes 



RmvA'^Jnl) ~ exp( yCT'S'g 



dS. 



6i'i, dlnv V 



(24) 



Both formulae have been shown to give reasonable agree- 
ment with numerical simulations of non- Gaussian cos- 
mologies [42, 109, 110] (but note that [12, 111] have 
reached somewhat different conclusions). Expanding 8i, 
at the first order in /^^ shows that these two theoretical 
expectations differ in the coefficient of the vaS-i term. 
Therefore, it is interesting to consider also the approxi- 
mation 



^(f^) /nl) ~ exp 



6 



'2^ ^ 6 dlni/ 



(25) 



which is designed to match better the Edgeworth ex- 
pansion of [107] when the peak height \s v ^ 1 [49]. 
When the primordial trispectrum is large (i.e. when 
5nl ~ 10^)1 terms involving the kurtosis must be in- 
cluded [49, 106, 107, 112]. In this case, it is also impor- 
tant to take into account a possible renormalization of the 
fiuctuation amplitude, ag — erg + 5a^ (Eq.l2), to which 
the high mass tail of the mass function is exponentially 
sensitive [49]. 

Figure 4 shows the effect of primordial NG of the lo- 
cal /jl^L type on the mass function of SO halos identified 
with a redshift-dependent overdensity threshold Avir(2) 
(motivated by spherical collapse in a ACDM cosmology 
[113]). Overall, the approximation Eq.(25) agrees better 
with the measurements than Eq.(24), which somewhat 
overestimates the data for f]^f^ ~ 100, and than Eq.(22), 
which is not always positive definite for f]^f^ = —100. 
However, as can be seen in Fig. 5, while the agreement 
with the data is reasonable for the SO halos, the theory 
strongly overestimates the effect in the mass function of 
EOF halos. Reference [110], who use a FOE algorithm 
with b ~ 0.2, makes the replacement 6c — ^ 6cy/q with 
q ~ 0.75 to match the Gaussian and non-Gaussian mass 
functions. A physical motivation of this replacement is 
provided by [114, 115], who demonstrate that the diffu- 
sive nature of the collapse barrier introduces a similar 
factor q = {1 + Db)~^ regardless the initial conditions. 
However, the value of the diffusion constant Db was ac- 
tually measured from simulations that use a SO finder 
with constant Avh- = 200 [116]. In the excursion set 
approach of [117], the value of q is obtained by normaliz- 
ing the Gaussian mass function to simulation (i.e. it has 
nothing to do with the collapse dynamics) and, there- 
fore, depends on the halo finder. Figure 5 demonstrates 
that this is also the case for the non-Gaussian correction 
i?(j^, /^l): choosing q ~ 0.75 as advocated in [110] gives 
good agreement for EOF halos, but strongly underesti- 
mates the effect for SO halos, for which the best-fit q is 
close to unity. As we will see below, the strength of the 
non-Gaussian bias may also be sensitive to the choice of 
halo finder. 
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More sophisticated formulations based on extended 
Press-Schechter (EPS) theory and/or modifications of 
the collapse criterion look promising since they can rea- 
sonably reproduce both the Gaussian halo counts and 
the dependence on [115, 118, 119]. The probability 
of first upcrossing can, in principle, be derived for any 
non-Gaussian density field and any choice of smoothing 
filter [120, 121]. For a general filter, the non-Markovian 
dynamics generates additional terms in the non-Gaussian 
correction to the mass function that arise from 3-point 
correlators of the smoothed density at different mass 
scales [115]. However, large error bars still make difficult 
to test for the presence of such sub-leading terms. For 
generic moving barriers B{a) such as those appearing in 
models of triaxial collapse [122, 123], the leading contri- 
bution to the non- Gaussian correction approximately is 
[118] 

Rii'jD^l + l'^S.H,^^) (26) 

where H^{v) = ~ ?>i> and the last equality assumes 
I' > 1. For SO halos, Eq.(26) with q ~ 0.7 does not fit to 
the measured correction Rlv, f^l) better than Eq.(25). 
However, the ellipsoidal collapse barrier with g ~ 0.7 
matches better the Gaussian mass function for moderate 
peak height <2 [119]. 

Parameterizations of the fractional correction based on 
N-body simulations have also been considered. While [43] 
considers a fourth-order polynomial fit to account for val- 
ues of large as 750, [12] exploits the fact that, for 
sufficiently small /^lj there is a one-to-one mapping be- 
tween halos in Gaussian and non-Gaussian cosmologies. 
In both cases, the fitting functions are consistent with 
the simulations at the few percent level. 

2. Clusters abundance 

Rich clusters of galaxies trace the rare, high-density 
peaks in the initial conditions and thus offer the best 
probe of the high mass tail of the multiplicity function. 
To infer the cluster mass function, the X-ray and mil- 
limeter windows are better suited than the optical-wave 
range because selection effects can be understood better 
(see, however, [124]). 

Following early theoretical [97, 98, 125-127] and nu- 
merical [36, 128-130] work on the effect of non-Gaussian 
initial conditions on the multiplicity function of cosmic 
structures, the abundance of clusters and X-ray counts 
in non-Gaussian cosmologies has received much atten- 
tion in the literature. At fixed normalization of the ob- 
served abundance of local clusters, the proto-clusters as- 
sociated with high redshift (2 < z < 4) Lya emitters are 
much more likely to develop in strongly non-Gaussian 
models than in the Gaussian paradigm [40, 111, 131]. 
Considering the redshift evolution of cluster abundances 




FIG. 5: Fractional deviation from the Gaussian mass func- 
tion as a function of the peak height v — Sc/cr. Different 
symbols refer to different redshifts as in Fig. 4. The curves 
are the theoretical prediction Eq. (25) at 2 = with q = 1 
(solid) and q — 0.75 (dotted). In the top panel, halos were 
identified using a spherical overdensity (SO) finder with a 
redshift-dependent overdensity threshold Avir{z) whereas, in 
the bottom panel, a Friends-of-Friends (FOF) finding algo- 
rithm with linking length 6 — 0.2 was used. 

thus can break the degeneracy between the initial den- 
sity PDF and the background cosmology. Simple exten- 
sions of the Press-Schechter formalism similar to those 
considered above have been shown to capture reasonably 
well the cluster mass function over a wide range of red- 
shift for various non-Gaussian scenarios [132]. Scaling 
relations between the cluster mass. X-ray temperature 
and Compton y-parameter calibrated using theory, ob- 
servations and detailed simulations of cluster formation 
[133, 134], can be exploited to predict the observed dis- 
tribution functions of X-ray and SZ signals and assess 
the capability of cluster surveys to test the nature of the 
initial conditions [135-141]. 

An important limitation of this method is that the 
impact of realistic models of primordial non-Gaussianity 
on cluster abundances is small compared to systematics 
in current and upcoming surveys [142-144]. Given the 
current uncertainties in the redshift evolution of clusters 
(one commonly assumes that clusters are observed at the 
epoch they collapse [143]), the selection effects in the cal- 
ibration of X-ray and SZ fluxes with halo mass, the free- 
dom in the definition of the halo mass, the degeneracy 
with the normalization amplitude erg (for positive /^li 
the mass function is more enhanced at the high mass 
end, and this is similar to an increase in the amplitude 



10 



of fluctuations ag [145]) and the low number statistics, 
the prospects of using the cluster mass function only to 
place competitive limits on with the current data 
are small. A two-fold improvement in cluster mass cali- 
bration is required to provide constraints comparable to 
CMB measurements [144]. 



the Gaussian model with identical mass power spectrum, 
subsequent work focused mostly on abundances (§IVB) 
or higher order statistics such as the bispcctrum (§IVD). 
It is only recently that [12] have demonstrated the strong 
scale-dependent bias arising in non-Gaussian models of 
the local f^l type. 



3. Voids abundance 



1. The non-Gaussian bias 



The frequency of cosmic voids offers a probe of the low 
density tail of the initial PDF [146]. The Press-Schechter 
formalism can also be applied to ascertain the sensitiv- 
ity of the voids abundance to non-Gaussian initial con- 
ditions. Voids are defined as regions of mass M whose 
density is less than some critical value ^„ < or, alterna- 
tively, whose three eigenvalues of the tidal tensor [147] lie 
below some critical value < [67, 119, 146, 148]. An 
important aspect in the calculation of the mass function 
of voids is the over-counting of voids located inside col- 
lapsing regions. This voids-in-clouds problem, as identi- 
fied by [149]), can be solved within the excursion set the- 
ory by studying a two barriers problem: 5c for halos and 
(5„ for voids. Including this eff'ect reduces the frequency 
of the smallest voids [119]. Neglecting this complication 
notwithstanding, the differential number density of voids 
of radius R is [146, 148] 



dn 
dR 



9 
2^ 



2 i?4 ^ ' dlnM 



l~aSsH3{\p,\) 



(27) 

where Vy = S^/cfm- While a positive /^l produces more 
massive halos, it generates fewer large voids [119, 146]. 
Hence, the effect is qualitatively different from a sim- 
ple rescaling of the normalization amplitude erg. A joint 
analysis of both abundances of clusters and cosmic voids 
might thus provide interesting constraints on the shape 
of the primordial 3-point function. There are, however, 
several caveats to this method, including the fact that 
there is no unique way to define voids [146]. Clearly, 
voids identification algorithms will have to be tested on 
numerical simulations [150] before a robust method can 
be applied to real data. 



C. Galaxy 2-point correlation 

Before [91] showed that, in Gaussian cosmologies, cor- 
relations of galaxies and clusters can be amplified rela- 
tive to the mass distribution, it was argued that primeval 
fluctuations have a non-Gaussian spectrum [151, 152] to 
explain the observed strong correlation of Abell clusters 
[153, 154]. Along these lines, [155] pointed out that pri- 
mordial non-Gaussianity could signiflcantly increase the 
amplitude of the two-point correlation of galaxies and 
clusters on large scales. However, except from [156] 
who showed that correlations of high density peaks in 
non-Gaussian models are signiflcantly stronger than in 



In the original derivation of [12], the Laplacian is ap- 
plied to the local mapping ^ = (j) + f^{^(j)^ in order to 
show that, upon substitution of the Poisson equation, 
the overdensity in the neighborhood of density peaks is 
spatially modulated by a factor proportional to the local 
value of (p. Taking into account the coherent motions in- 
duced by gravitational instabilities, the scale-dependent 
bias correction reads 



A6«(fc, f^l) = 3f^l [h{M) - 1] 5c{0) 



k^T{k)D{z) 



, (28) 



where bi{M) is the linear, Eulerian bias of halos of mass 
M . This effect can be understood intuitively in terms of a 
local rescaling of the small-scale amplitude of matter fluc- 
tuations or, equivalently, a local rescaling of the critical 
density threshold [12, 157]. The original result missed out 
a multiplicative factor of T(k)~^ which was introduced 
subsequently by [158] upon a derivation of Eq. (28) in 
the limit of high density peaks. The peak-background 
split approach [104, 159, 160] promoted by [157] shows 
that the scale-dependent bias applies to any tracer of the 
matter density held whose (Gaussian) multiplicity func- 
tion depends on the local mass density only. In this ap- 
proach, the Gaussian piece of the potential is decomposed 
into short- and long- wavelength modes, 4> = (f)i + (f)s. The 
short-wavelength piece of the density fleld is then given 
by the convolution 

5s = Mi.^s = M^cPs (1 + VNl<Pi)+f^"lM*<f>l , (29) 

where M. is the transfer function Eq.(7). Ignoring the 

white-noise t(Tm, this provides an intuitive explanation 
of the efl'ect in terms of a local rescaling of the small-scale 
amplitude of matter fluctuations. 



(30) 



Assuming the mass function depends only on the peak 
height v = 5c/ as, the long- wavelength part of the halo 
overdensity becomes [157] (see also [12, 72, 161] 



1 



/loc 
NL' 



n(zv) \^(7s(l + 2/j 

-l(5ii^)+2f^lcPi{^)) 
cr,, \ / 




1 



(31) 



dv 



Upon a Fourier transformation and using the fact that, in 
the Gaussian case, 5f{]s.) = 6L<5i(k) with the Lagrangian 
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FIG. 6: Non-Gaussian bias correction Eq.(33) for the local, 
equilateral and folded primordial bispectrum as a function of 
wavenumber . Results are shown at z = for a smoothing 
radius R = 5 h~^Mpc and a nonlinear parameter /^l ~ 1- 
The dotted line represent an analytic approximation, Eqs (35) 
- (37), which is valid at large scales A; ^ 1. Note that, 
while the magnitude of Ab^ik, /nl) does not change with R, 
Abnik, /nl) ^■iid Abnik, /nl,) strongly depends on the smooth- 
ing radius (see text). 



bias 6l — —cr^^dlnn/di', we recover the non-Gaussian 
bias correction Eq.(28) provided that the tracers move 
coherently with the dark matter, i.e. &l = bi{M) — 1 
[162]. As emphasized in [12], the scale-dependence arises 
from the fact that the non-Gaussian curvature pertur- 
bation $(x) is related to the density through the Pois- 
son equation (6) (so that (5;(k) = M{k)(f>i{'k)). There is 
no such effect in the (local) model, Eq.(lO), nor in 
texture-seeded cosmologies [163] for instance. 

The derivation of [158], based on the clustering of re- 
gions of the smoothed density field 5m above threshold 
dc{z), is formally valid for high density peaks only. How- 
ever, it is general enough to apply to any shape of pri- 
mordial bispectrum. The 2-point correlation function of 
that level excursion set, which was first derived by [125], 
can be expressed in the high threshold limit (z^ ^ 1) as 



oo n — 1 



^>A^) = -1 + exp<^ — — 

I j\[n- jy. 



(32) 



x;i, ■ ■ • ,Xi, 
j times 



X2, • • 

(" - j) 



,x;2 
times 



,z = 



where r = xi — X2. For most non-Gaussian models in 
which the primordial 3-point function is the dominant 
correction, this expansion can be truncated at the third 
order and Fourier transformed to yield the non-Gaussian 
correction AP>^(fc) to the power spectrum. Assum- 
ing a small level of primordial NG, we can also write 



AF>^(fc) « 2bLAb^PR{k) where &l 
tually obtain 



ly'^/Sc, and even- 



A6„(fc,/^L)s64fc)^(fc,/^L) 



Xfl(fc,0) 



|-^(fc,/NL). 

(33) 

The dependence on the shape of the 3-point function is 
encoded in the function J^(fc, /^l) [158, 164], 



dki klMniki^O) 



(34) 



/^^j »^ / I- ,^^B^{ki,^/a,k) 
X / d^A^flXVa, 0)- 



P*(fc) 



where = P + kf + 2^kki. Note that, for f^l < 
0, this first order approximation always breaks down at 
sufficiently small k because AP>^(fc) < 0. 

Figure 6 shows the non-Gaussian halo bias Eq.(33) in- 
duced by the local, equilateral and folded bispectrum 
[164] . In the local and folded non-Gaussianity, the de- 
viation is negligible at k — 0.1 /iMpc"^, but increases 
rapidly with decreasing wavenumber. Still, the large 
scale correction is much smaller for the folded template, 
and nearly absent for the equilateral type, which make 
them much more difficult to detect with galaxy surveys 
[164]. To get insights into the behavior of lS.h,^{k^ /nl) 
large scales, let us identify the dominant contribution to 
T{k,f^-^) in the limit A: < 1. Setting MR{y/a,0) w 
^Afl{kl,0) and expanding P^{y/a) at second order in 
k/ki, we find after some algebra 



-^(^7 /nl) 
•^(^' /nl) 



•^(^7 /nl) 



r\oc 

/nl 



3S 



2(n, - 4)\ 2(4-„,) 



.(M) 



k^r 



+ -(n,-4)I]J-2)fc2 

3 rfol V. f'^s - 4\ 



'R 
-2 

R ' 



assuming no running scalar index, i.e. dng/dlnk 
The auxiliary function S^(?t.) is defined as 



1 

^ dfcfc(2+«)Xfl(fc,0)2p4fc) 



(35) 
(36) 

(37) 
= 0. 

(38) 



Hence, we have 2^(0) = '^r- '^^^ be seen in Fig. 6, 
these approximations capture relatively well the large 
scale non- Gaussian bias correction induced by the equi- 
lateral and folded type of non-Gaussianity. For a nearly 
scale-invariant spectrum « 1, the effect scales as 
Abu oc k and A6k cx const., respectively. Another 
important feature of the equilateral and folded non- 
Gaussian bias is the dependence on the mass scale M 
through the multiplicative factor cr^^. Indeed, choosing 
R = 1 h~^Mpc instead of i? = 5 h~^Mpc as done in 
Fig. 6 would suppress the effect by a factor of ^ 3. In 



the high peak limit, cr^ 



bh/^ciz) which cancels out 



the dependence on redshift but enhances the sensitivity 



12 



to the halo bias, i.e. Ab^ oc 6^ for the equilateral and 
folded shapes whereas Ab^ (x bi, in the local model. 

At this point, it is appropriate to mention a few caveats 
to these calculations. Firstly, Eq. (28) assumes that 
the tracers form after a spherical collapse, which may 
be a good approximation for the massive halos only. If 
one instead considers the ellipsoidal collapse dynamics, in 
which the evolution of a perturbation depends upon the 
three eigenvalues of the initial tidal shear, (5c (0) should 
be replaced by its ellipsoidal counterparts (5oc(0) which 
is always larger than the spherical value [122]. In this 
model, the scale-dependent bias Afe^ is thus enhanced by 
a factor (5oc(0)/(5c(0) [12, 161]. Secondly, Eq. (28) assumes 
that the biasing of the surveyed objects is described by 
the peak height u only or, equivalently, the hosting halo 
mass M . However, this may not be true for quasars 
whose activity may be triggered by merger of halos [165, 
166]. Reference [157] used the EPS formalism to estimate 



the bias correction Ah-n 
Ab 



merge: 



induced by mergers. 



(39) 



so the factor bi {M) — 1 should be replaced by bi {M) — 1 — 
« bi{M) - 1.6. The validity of this result should be 
evaluated with cosmological simulations of quasars for- 
mation. In this respect, semi-analytic models of galaxy 
formation suggest that merger-triggered objects such as 
quasars do not cluster much differently than other trac- 
ers of the same mass [167]. However, this does not mean 
that the same should hold for the non-Gaussian scale de- 
pendent bias. Still, since the recent merger model is an 
extreme case it seems likely that the actual bias correc- 
tion is < A6merger < . Thirdly, the scale-dependent 
bias has been derived using the Newtonian approxima- 
tion to the Poisson equation, so one may wonder whether 
general relativistic (GR) corrections to M.R{k)^^ may 
suppress the effect on scales comparable to the Hubble 
radius. Reference [168] showed how large scale primor- 
dial NG induced by GR corrections propagates onto small 
scales once cosmological perturbations reenter the Hub- 
ble radius in the matter dominated era. This effect gen- 
erates a scale-dependent bias comparable, albeit of op- 
posite sign to that induced by local NG [164]. More re- 
cently, [169] argues that there are no GR corrections to 
the non-Gaussian bias and that the scaling Ab^, (x k^^ 
applies down to smallish wavenumbers. 

We can also ask ourselves whether higher-order terms 
in the series expansion (32) furnish corrections to the 
non-Gaussian bias similar to Eq.(28). The quadratic cou- 
pling /nl'/'^ induces a second order correction to the halo 



power spectrum which reads [49] 



APh(fc) = -{fh°l 



bi{M) - 1 
X MR{k,{))P^{k) . 



'5l{z)s'i\M) 



(40) 



Its magnitude relative to the term linear in /^^i Eq.(28), 
is approximately 0.03 at redshift z = 1.8 and for a halo 
mass M = 10^3 Mq//i. Although its contribution be- 
comes increasingly important at higher redshift, it is 
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FIG. 7: Non-Gaussian bias correction (filled symbols) for ha- 
los of mass M > 2 X 10^^ Mq/H extracted at z = 0.5 from 
simulations of the local /^l model. The solid curve repre- 
sents the theoretical model Eq. (44). The dashed, dotted- 
dashed-dotted and dotted curves show the scale-dependent 
bias A6k, the scale-independent oflFset A6i and the contribu- 
tion from the matter power spectrum bi{M)l3s that arise at 
first order in fj^l (see Fig.2). The error bars indicate the 
scatter among 8 realizations of 1024'' simulations with box 
size L = 1600 /i"^Mpc. 



fairly small for realistic values of /^l • In local NG model, 
the power spectrum of biased tracers of the density field 
can also be obtained from a local Taylor series in the 
evolved (Eulerian) density contrast 6 and the Gaussian 
part (f> of the initial (Lagrangian) curvature perturbation 
[48, 72]. Using this approach, it can be shown that the 
halo power spectrum arising from the first order terms of 
the local bias expansion can be cast into the form [48] 

Ph(fc) = HM) + f^°lb4k)] 'Pnik) (41) 

Hence, we also obtain a second order term propor- 
tional to {f:^l)^M^^PR{k) = {f^l)^P^{k) which, how- 
ever, contributes only at very small wavenumber k < 
0.001 h-^Mpc. All this suggests that Eq. (28) is the 
dominant contribution to the non-Gaussian bias in the 
wavenumber range 0.001 ^ k < 0.1 /iMpc~^. 

Finally, a non-Gaussian, scale-dependent bias correc- 
tion can also arise in the local, deterministic bias ansatz 

(5h(x) = 6i(5(x)+62<5(x)V2H [170] if the initial density 

field is non-Gaussian. Here, b^ is the A^th-order bias pa- 
rameters (here again, the first-order bias is 6i = 1 + &l)- 
In this approach, the correction is induced by the correla- 
tion foi62(<5(xi)<5^(x2)) between the linear and quadratic 
term in the galaxy biasing relation (which is in fact a 
collapsed or squeezed 3-point function) and thus reads 
[68, 69] 



A6.(fc,/;°£) = 2f^°lb2aj,MR{k,0)-^ 



(42) 
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Even though 62 (t|. « &l<5c in the high-threshold hmit 
1/ ^ 1, b2<jji behaves very differently than 6l(5c for mod- 
erate peak height because 62 is proportional to the sec- 
ond derivative of the mass function n^v). So far however, 
Eq.(28) appears to describe reasonably well the numeri- 
cal results for a wide range of halo bias. 



2. Comparison with simulations 

In order to fully exploit the potential of forthcoming 
large-scale surveys, a number of studies have tested the 
theoretical prediction against the outcome of large nu- 
merical simulations [12, 42-44, 72, 110]. 

At the lowest order, there are two additional albeit 
relatively smaller corrections to the Gaussian bias which 
arise from the dependence of both the halo number den- 
sity n{M,z) and the matter power spectrum Ps{k,z) 
on primordial NG [42]. Firstly, assuming the peak- 
background split holds, the change in the mean num- 
ber density of halos induces a scale-independent offset 
which we denote A6i(/^l)- terms of the non-Gaussian 
fractional correction Riv, /nl) to the mass function, this 



contribution is 



A6i(/, 



loc\ 
NL j 



1 d 



di 



In 



(43) 



It is worth noticing that Abi{fj^j^) has a sign opposite 
to that of because the bias decreases when the 

mass function goes up. In practice, the approximation 
Eq. (25), which matches well the SO data for ^ 4, can 
be used for moderate values of the peak height. For FOF 
halos with linking length b — 0.2, one should make the 
replacement S^. Scy/q with q « 0.75 in the calculation 
of the scale- independent offset. It is sensible to evaluate 
A6i(/^l) ^t ^ mass scale (M) equal to the average halo 
mass of the sample. Secondly, we also need to account 
for the change in the matter power spectrum (see Fig. 2 
in §111). Summarizing, local non-Gaussianity adds a cor- 
rection Ab{k, /nl) to the bias b{k) of dark matter halos 
that reads [42] 

A6(fc, f^°l) = A6«(fc, f^°l) + A6i(/;°£) + 6i(M)/35(fc, fh°l) 

(44) 

at first order in f^£. As can be seen in Fig. 7, the inclusion 
of these extra terms substantially improves the compari- 
son between the theory and the simulations. Considering 
only the scale-dependent shift Ab^ leads to an appar- 
ent suppression of the effect in simulations relative to 
the theory. Including the scale-independent offset Abi 
considerably improves the agreement at wavenumbers 
k < 0.05 /iMpc~^. Finally, adding the scale-dependent 
term bi{M)f3s further adjusts the match at small scale 
k > 0.05 hMpc~^ by making the non-Gaussian bias shift 
less negative. Along these lines, [72] find that the inclu- 
sion of Abi to the bias also improves the agreement with 
measurements of Ab{k, f^£) obtained for FOF halos, and 
show that taking into account second- and higher-order 
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FIG. 8: Halo-halo (solid curve) and halo-matter (dashed 
curve) power spectra Ph{k) and Pi^s{k) measured in simu- 
lations of the Gaussian and /nl = ±100 models for ha- 
los of mass M > 2 X lO" Mq/H at redshift z = 1. The 
error bars represent the scatter among 8 realizations. For 
/nl ~ —100, the cross-power spectrum is negative on scales 



k < 0.005 /iMpc" 
prediction. 



in good agreement with the theoretical 



corrections could extend the validity of the theory up to 
scales fc ~ 0.1 - 0.3 /iMpc"^ 

The non-Gaussian bias correction can be measured in 
the cross- and auto-power spectrum of dark matter ha- 
los, Phsik) and Ph{k). To compute these quantities, 
dark matter particles and halo centers are interpolated 
onto a regular cubical mesh. The resulting dark matter 
and halo fluctuation fields, (5m (k) and (5h(k), are then 
Fourier transformed to yield the matter-matter, halo- 
matter and halo-halo power spectra Ps{k), Phs{k) and 
Ph{k), respectively. Ph{k) is then corrected for the shot 
noise, which is assumed to be l/fih if dark matter halos 
are a Poisson sampling of some continuous field. This 
discreteness correction is negligible for Ps{k) due to the 
large number of dark matter particles. On linear scales 
(fc < 0.01 /iMpc-i), the halo bias 6(k) = (5h(k)/5m(k) 
approaches the constant value bi{M) which needs to be 
measured accurately as it controls the strength of the 
scale-dependent bias correction A6„. In this respect, the 
ratio Phs{k)/Ps{k) is a better proxy for the halo bias 
since it is less sensitive to shot-noise. 

Auto- and cross-power analyses may not agree with 
each other if the halos and dark matter do not trace 
each other on scale k < 0.01 hMpc~^ where the non- 
Gaussian bias is large, i.e. if there is stochasticity. Fig. 8 
shows Phs{k) and Ph{k) averaged over 8 realizations of 
the models with f^^l = 0,±100 [42]. The same Gaus- 
sian random seed field (f> was used in each set of runs 
so as to minimize the sampling variance. Measurements 



14 



of the non-Gaussian bias correction obtained with the 
halo-halo or the halo-matter power spectrum are in a 
good agreement with each other, indicating that non- 
Gaussianity does not induce stochasticity and the pre- 
dicted scaling Eq.(28) applies equally well for the auto- 
and cross-power spectrum. However, while a number of 
numerical studies of the model have confirmed the 
scaling A/;,-(/c, oc Mii{k)~^ and the redshift depen- 
dence oc D{z)~^ [12, 42, 43, 110], the exact amplitude of 
the non- Gaussian bias correction remains somewhat de- 
batable. Reference [42] who use SO halos and [72] who 
use FOF halos find satisfactory agreement with the the- 
ory once the scale-independent offset A6i is included. By 
contrast, [43], who use the same FOF halos as [72], argue 
that the scale-dependent piece Ab^ requires, among oth- 
ers, a multiplicative correction of the form (1 — 
with /3i ^ 4 X 10"'* > 0. Similarly, [110] who also use 
FOF halos find that the theory is a good fit to the sim- 
ulations only upon replacing 6l by qbh in Eq.(33), with 
q 0.75. Part of the discrepancy may be probably due 
to the fact that the last two references do not include 
A6i, which leads to an apparent suppression of the effect 
(see Fig. 7). Another possible source of discrepancy may 
be choice of the halo finder which, as seen in Fig. 5, has 
an impact on the strength of the non-Gaussian correction 
to the mass function. This possibility is investigated in 
Figure 9, which shows the non-Gaussian bias correction 
obtained with FOF halos. For this low biased sample, 
the scale- independent correction is |A6i| < 0.003 and 
can thus be neglected. The best-fit values of sig- 
nificantly below the input values of ±100, in agreement 
with the findings of [43, 110] (note, however, that this 
suppression is more consistent with Sc being rescaled by 
^/q6c « 0.865c and &l being unchanged). This indicates 
that the choice of halo finder may also affect the mag- 
nitude of the scale-dependent non-Gaussian bias. Dis- 
crepancies have also been observed between the theoret- 
ical and measured non-Gaussian bias corrections in non- 
Gaussian models of the local cubic-order coupling g]^l^4>^ 
[49]. Understanding all these results clearly requires a 
better modeling of halo clustering. 



3. Redshift- space distortions 

Peculiar velocities generate systematic differences be- 
tween the spatial distribution of data in real and redshift 
space. These redshift-space distortions must be properly 
taken into account in order to extract from redshift 
surveys. On the linear scales of interest, the redshift- 
space power spectrum of biased tracers reads as [171, 172] 

P^{k,fi)= [blPs{k)+2bif^^Ps0{k)+ffi''Pe{k)\ , (45) 

where Pse and Pg are the density-velocity and velocity 
divergence power spectra, n is the cosine of the angle 
between the wavemode k and the line of sight and / is 
the logarithmic derivative of the growth factor. For Pg, 



the 1-loop correction due to primordial NG is identical 
to Eq.(17) provided F2(ki,k2) is replaced by the kernel 

G2(ki,k2) = 3/7+/i(fci/fc2 + fc2/fci)/2 + 4^V7 describing 
the 2nd order evolution of the velocity divergence [63]. 
For Pse, this correction is 



A pNG 



(k) 



f^2(q,k-q) + G2(q,k-q) 



(27r)3 
X Bo(-k,q,k-q) . 



(46) 



Again, causality implies that G2(ki,k2) vanishes in the 
limit ki = — k2. For unbiased tracers with bi = 1, the 
linear Kaiser relation is thus recovered at large scales 
k < 0.01 /iMpc~* (this is consistent with the analysis of 
[173]). For biased tracers, we still expect the Kaiser for- 
mula to be valid, but the distortion parameter /3 should 
now be equal to /3 = f/{bi + Ab^), where Abn{k, /^l) is 
the scale-dependent bias induced by the primordial non- 
Gaussianity. 



4- Mitigating cosmic variance and shot-noise 

Because of the finite number of large scale wavemodes 
accessible to a survey, any large scale measurement of 
the power spectrum is limited by the cosmic (or sam- 
pling) variance caused by the random nature of the wave- 
modes. For discrete tracers such as galaxies, the shot 
noise is another source of error. Restricting ourselves 
to weak primordial NG, the relative error on the power 
spectrum P is ap/P ss 1/Vn{1 + al/P), where N is 
the number of independent modes measured and ct^ is 
the shot-noise [174]. Under the standard assumption of 
Poisson sampling, cr^ equals the inverse of the number 
density 1/n and causes a scale- independent enhancement 
of the power spectrum. The extent to which one can im- 
prove the observational limits on the nonlinear parame- 
ters ''^ill strongly depend on our ability to minimize 
the impact of these two sources of errors. By comparing 
differently biased tracers of the same surveyed volume 
[175, 176] and suitably weighting galaxies (by the mass 
of their host halo for instance) [177, 178], it should be 
possible to circumvent these problems and considerably 
improve the detection level. 

Figure 9 illustrates how the impact of sampling vari- 
ance on the measurement of can be mitigated. 
Namely, the data points show the result of taking the ra- 
tio Ph(fc, fj^l)/Ps{k, f^£) for each set of runs with same 
Gaussian random seed field (f> before averaging over the 
realizations. This procedure is equivalent to the multi- 
tracers method advocated by [175]. Here, Ps can be 
thought as mimicking the power spectrum of a nearly 
unbiased tracer of the mass density field with high num- 
ber density. Although, in practical applications, using 
the dark matter field works better [179], in real data Ps 
should be replaced by a tracer of the same surveyed vol- 
ume different than the one used to compute Ph- Figure 9 
also shows that, upon taking out most of the cosmic vari- 
ance, there is some residual noise caused by the discrete 
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FIG. 9: Fractional correction to the Gaussian halo bias in the 
/nl ~ ±100 and Gaussian models. In contrast to Fig. 8, halos 
were identified with a FOF finder of linking length b — 0.2. 
Only the wavemodes to the left of the vertical line were used 
to fit Abfi{k, f^£). For this low biased sample, the scale- 
independent correction is |A6i| < 0.003 and can thus be ig- 
nored. The best-fit value of /J^l and the corresponding Icr 
error is quoted for each model (Figure taken from [179]). 

nature of the dark matter halos. As shown recently [178] 
however, weighting the halos according to their mass can 
dramatically reduce the shot noise relative to the Pois- 
son expectation, at least when compared against the dark 
matter. Applying such a weighting may thus significantly 
improve the error on the nonlinear parameter /^li but 
this should be explored in realistic simulations of galax- 
ies, especially because the halo mass M may not be easily 
measurable from real data [179] . This approach undoubt- 
edly deserves further attention as it has the potential to 
substantially improve the extraction of the primordial 
non-Gaussian signal from galaxy surveys. 

To conclude this Section, it is worth noting that, while 
the PDF of power values P(k) has little discriminatory 
power (for large surveyed volume, it converges towards 
the Rayleigh distribution as a consequence of the cen- 
tral limit theorem) [180], the covariance of power spec- 
trum measurements (which is sensitive to the selection 
function, but also to correlations among the phase of the 
Fourier modes) may provide quantitative limits on cer- 
tain type of non-Gaussian models [174, 181]. 



D. Galaxy bispectrum and higher order statistics 

Higher statistics of biased tracers, such as the galaxy 
bispectrum, are of great interest as they are much more 
sensitive to the shape of the primordial 3-point function 
than the power spectrum [11, 44, 69, 182, 183]. There- 
fore, they could break some of the degeneracies affecting 
the non-Gaussian halo bias (For example, the leading 
order scale-dependent correction to the Gaussian bias in- 



duced by the local quadratic and cubic coupling are fully 
degenerated [49]). 

1. Normalized cumulants of the galaxy distribution 

The skewness of the galaxy count probability distribu- 
tion function could provide constraints on the amount 
of non-Gaussianity in the initial conditions. As dis- 
cussed in §111 however, it is difficult to disentangle the 
primordial and gravitational causes of skewness in low 
redshift data unless the initial density field is strongly 
non-Gaussian. The first analyzes of galaxy catalogs in 
terms of count-in-cells densities all reached the conclu- 
sion that the skewness (and higher-order moments) of the 
observed galaxy count PDF is consistent with the value 
predicted by gravitational instability of initially Gaus- 
sian fluctuations [51, 58, 61, 184-186]. Back then how- 
ever, most of the galaxy samples available were not large 
enough to accurately determine the Sj at large scales 
[187]. Despite the two orders of magnitude increase in 
surveyed volume, these measurements are still sensitive 
to cosmic variance, i.e. to the presence of massive super- 
clusters or large voids. Nevertheless, the best estimates 
of the first normalized cumulants Sj of the galaxy PDF 
strongly suggest that high order galaxy correlation func- 
tions indeed follow the hierarchical scaling predicted by 
the gravitational clustering of Gaussian ICs [188]. There 
is no evidence for strong non-Gaussianity in the initial 
density field as might by seeded by cosmic strings or tex- 
tures [189]. 

The genus statistics of constant density surfaces 
through the galaxy distribution measures the relative 
abundance of low and high density regions as a func- 
tion of the smoothing scale R and, therefore, could 
also be used as a diagnostic tool for primordial non- 
Gaussianity. For a Gaussian random field, the genus 
curve (i.e. the genus number as a function of the density 
contrast) is symmetric about 5/? = regardless the value 
of R. Primordial NG and nonlinear gravitational evolu- 
tion can disrupt this symmetry [190]. The effect of non- 
Gaussian ICs on the topology of the galaxy distribution 
has been explored in a number of papers [36, 191-194]. 
For large values of R and realistic amount of primordial 
NG, the genus statistics can also be expanded in a se- 
ries whose coefficients are the normalized cumulants Sj 
of the smoothed galaxy density field. In other words, the 
genus statistics essentially provides another measure of 
the (large scale) cumulants. So far, measurements from 
galaxy data are broadly consistent with Gaussian initial 
conditions [195, 196]. 

2. Galaxy bispectrum 

Most of the scale-dependence of the primordial n-point 
functions is integrated out in the normalized cumulants, 
which makes them weakly sensitive to primordial NG. 
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However, while the eflect of non-Gaussian initial con- 
ditions, galaxy bias, gravitational instabilities etc. are 
strongly degenerated in the Sj, they imprint distinct sig- 
natures in the galaxy bispectrum i?g(ki, k2, ka), an ac- 
curate measurement of which could thus constrain the 
shape of the primordial 3-point function. 

In the original derivation of [182], the large scale (un- 
filtered) galaxy bispectrum in the model is given by 



Sg(ki,k2,k3) = &?So(ki,k2,k3) 



(47) 



+ 6?&2 Po(fci)Po(fc2) + (cyc.) 



26? F2(ki,k2)Po(A:i)Po(fc2) + (cyc.) 



Again, bi and &2 are the first- and second-order bias 
parameters that describe the galaxy biasing relation as- 
sumed local and deterministic [170]. The first term in the 
right-hand side is the primordial contribution which, for 
equilateral configurations and in the fj^f^ model, scales 
as Mnik, z)~^ like in the matter bispectrum, Eq.(18). 
The two last terms are the contribution from nonlinear 
bias and the trcc-lcvcl correction from gravitational in- 
stabilities, respectively. They have the smallest signal in 
squeezed configurations. 

As recognized by [69, 183], Eq.(47) misses an impor- 
tant term that may significantly enhance the sensitivity 
of the galaxy bispectrum to non-Gaussian initial condi- 
tions. This contribution is sourced by the trispectrum 
TR(ki, k2, ks, k4) of the smoothed mass density field, 



q) + (2 perms.) . (48) 



At large scale, this simplifies to the sum of the lin- 
early evolved primordial trispectrum To(ki, k2, ka, k4) 
and a coupling between the primordial bispectrum 
i3o(ki. k2, ks) (linear in /^l) a-^id the second order PT 
corrections (through the kernel F2(ki, k2)). In the case of 
local non-Gaussianity and for equilateral configurations, 
the first piece proportional to To scales as (/nl)^^^^ 
times the Gaussian tree-level prediction, with the same 
redshift dependence. Hence, it is similar to the second 
order correction {f^i^^ M.^ Pnik) that appears in the 
halo power spectrum (see Eq.41). The second piece lin- 
ear in generates a signal at large scales for essen- 
tially all triangle shapes in the local model as well as in 
the case of equilateral NG. This second contribution is 
maximized in the squeezed limit (where it is one order 
of magnitude larger than the result obtained by [182]) 
which helps disentangling it from the Gaussian terms. 
Note that a strong dependence on triangle shape is also 
present in other NG scenarios such as the model [63]. 

This newly derived contributions are claimed to lead 
to more than one order of magnitude improvement in 
certain limits [183], but it is not yet clear whether these 
gains can be fully realized with upcoming galaxy sur- 
veys. To accurately predict the constraints that could be 
achieved with future measurements of the galaxy bispec- 



trum, a comparison of these predictions with the halo bis- 
pectrum extracted from numerical simulations is highly 

desirable. To date, the only numerical study [44] has 
measured the halo bispectrum for some isosceles triangles 
(fci = fe). While the shape dependence is in reasonable 
agreement with the theory, the observed fc-dependence 
appears to depart from the predicted scaling. 



E. Intergalactic medium and the Lya forest 

Primordial non-Gaussianity also affects the intergalac- 
tic medium (IGM) as a positive enhances the forma- 
tion of high-mass halos at early times and, therefore, ac- 
celerates reionization [197-199]. At lower redshift, small 
box hydrodynamical simulations of the Lya forest indi- 
cate that non-Gaussian initial conditions could leave a 
detectable signature in the Lya flux PDF, power spec- 
trum and bispectrum [200]. However, while differences 
appear quite pronounced in the high transmissivity tail 
of the flux PDF (i.e. in underdense regions), the Lya 
ID flux power spectrum seems little affected. Given the 
small box size of these hydrodynamical simulations, it is 
worth exploring the effect in large N-body cosmological 
simulations using a semi-analytic modeling of the Lya 
forest [201], even though such an approach only provides 
a very crude approximation to the temperature-density 
diagram of the IGM in hydrodynamical simulations. Fig- 
ure 10 shows the imprint of local type NG on the Lya 
3D flux power spectrum (which is not affected by pro- 
jection effects) extracted at z = 2 from a series of large 
simulations. The Lya transmitted flux is calculated in 
the Gunn-Peterson approximation [202]. A clear signa- 
ture similar to the non-Gaussian halo bias can be seen. 
As expected, it is of opposite sign since the Lya forest is 
anti-biased relative to the mass density ficild (overdensi- 
ties are mapped onto relatively low flux transmission). 

To estimate the strength of the signal (see [201] for 
the details), one can assume that the (real space) optical 
depth r(x) to Lya absorption at comoving position x is 
approximately [203] 



t(x) =f l + (5g(x) 



(49) 



where 5g is the gas density, r ~ 1 is the optical depth at 
mean gas density and a ~ 1 — 2 is some parameter that 
depends on the exact thermal history of the low density 
IGM. The above relation holds for the moderate over- 
densities 5g < 10 that are responsible for most of the 
Lya absorption features. To relate the gas density to the 
smoothed linear density field, we could make the simple 
ansatz 5g = Sr [204]. In this linear approximation how- 
ever, the large scale bias bp of the Lya flux density fleld 
is much larger than that measured in detailed numerical 
simulations (e.g., &| ~ 0.017 at z = 3 [205]). There- 
fore, one may want to consider the lognormal mapping 
[206, 207] 



1 + Sg = exp((5_R - cr|/2) 



(50) 



17 



1.4 
1.3 
1.2 
1.1 

1 

0.9 
0.8 
0.7 
0.6 



fnl = -100, resolution = 256"2 
fnl = +100, resolution = 256"2 
1 



+ +++ 



XX "J "J"® ><■ i ******** * * * 



0.01 



0.1 



k (h/Mpc) 



FIG. 10: Ratio between the z = 2 Lya flux power spectrum 
extracted from simulations of Gaussian and non-Gaussian ini- 
tial conditions. The mean transmission is set to F = 0.8 and 
the power-law exponent a = 1.65 (see text). 



to better capture nonlinearities in the gas density field. 
Expanding exp((5ii;) at second order and noticing that, 
in the presence of weak non-Gaussianity, the joint PDF 
P(5j^(xi), 5fl(x2)) can generically be expanded into an 
Edgeworth series where the primordial 3-point function is 
the dominant correction, it is straightforward to compute 
the Lya 3D flux power spectrum for nonzero /^l- Upon 
a Fourier transformation, we arrive at 



PF(fc,/^L) 
PF(fc,0) 



l-4gF4-^ii(fc,z)-ij-(fc,/] 



(51) 



where is some auxiliary function of (r, a, a^)- This re- 
sult is valid for any model of primordial NG characterized 
by an initial bispectrum. In the model, the large- 
scale non-Gaussian Lya bias scales as Ab-p{k, f^^J « 
—2g-paj^MR{k,z)^^ (x k^'^T{k)~^ like the non-Gaussian 
halo bias. Assuming f — 0.7, = 1.8 and a = 
1.65 yields a mean flux F w 0.8 and a ratio Ppik = 
0.01,/^L)/PF(fc = 0.01,0) « 1 T 0.13 for f^l = ±100 
comparable in magnitude to that seen in Fig. 10. A 
detection of this effect, although challenging in particu- 
lar because of continuum uncertainties, could be feasible 
with future data sets. Summarizing, the Lya should pro- 
vide interesting information on the non-Gaussian signal 
over a range of scale and redshift not easily accessible to 
galaxy and CMB observations [200, 201]. 



parameter(s) with LSS is to use large scale, high-redshift 
observations [331. 



A. Existing constraints on primordial NG 

The non-Gaussian halo bias presently is the only LSS 
method that provides a robust limit on the magnitude 
of a primordial 3-point function of the local shape. It 
is a broadband effect that can be easily measured with 
photometric redshifts. The authors of [157] have applied 
Eq.(28) to constrain the value of using a compila- 
tion of large-scale clustering data. Their constraint arise 
mostly from the QSO sample at median redshift z — 1.8, 
which covers a large comoving volume and is highly bi- 
ased, bi — 2.7. They obtain 



-29 < fl,°{: < +69 



(52) 



at 95% confidence level. These limits are competitive 
with those from CMB measurements, — 10 < < +'''4 
[208]. It is straightforward to translate this 2-(t limit 
into a constraint on the cubic order coupling yj^^ since 
the non-Gaussian scale-dependent bias A5k(/c,5nl) 
the same functional form as A6„(fc, J^^) [49]. Assuming 
/nl = 0' oue obtains 



-3.5 X 10^ < g]^l < +8.2 x 10^ 



(53) 



These limits are comparable with those inferred from the 
analysis of CMB data. 

Measurements of the galaxy bispectrum in several red- 
shift catalogs have shown evidence for a configuration 
shape dependence in agreement with that predicted from 
gravitational instability, ruling out initial conditions 
at the 95% C.L. [209, 210]. Recent analyses of the 
SDSS LRGs catalogue indicate that the shape depen- 
dence of the reduced 3-point correlation Q3 ^ $3/(^2)^ is 
also consistent with Gaussian ICs [211], although a pri- 
mordial (hierarchical) non-Gaussian contribution in the 
range Q3 ~ 0.5 — 3 cannot be ruled out [212]. Other 
LSS probes of primordial non-Gaussianity, such as the 
abundance of massive clusters, are still too affected by 
systematics to furnish tight constraints on the shape and 
magnitude of a primordial 3-point function. Still, the 
observation of a handful of unexpectedly massive high- 
redshift clusters has been interpreted as evidence of a 
substantial degree of primordial NG [213-215]. 



B. Future prospects 



V. CURRENT LIMITS AND PROSPECTS 

As the importance of primordial non-Gaussianity rela- 
tive to the non-Gaussianity induced by gravitational clus- 
tering and galaxy bias increases towards high redshift, 
the optimal strategy to constrain the nonlinear coupling 



Improving the current limits will further constrain the 
physical mechanisms for the generation of cosmological 
perturbations. 

The non-Gaussian halo bias also leaves a signature in 
cross-correlation statistics of weak cosmic shear (galaxy- 
galaxy and galaxy-CMB) [216, 217] and in the integrated 
Sachs- Wolfe (ISW) effect [157, 161, 218]. Measurements 
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of the lensing bispectrum could also constrain a number 
of non-Gaussian models [219]. However, galaxy cluster- 
ing will undoubtedly offer the most promising LSS di- 
agnostic of primordial non-Gaussianity. The detectabil- 
ity of a local primordial bispectrum has been assessed 
in a series of papers. It is expected that future all-sky 
galaxy surveys will achieve constraints of the order of 
^/nl ~ 1 assuming all systematics are reasonably under 
control [48, 107, 157, 161, 218, 220- 222]. Realistic mod- 
els of cubic type non-Gaussianity [49], modifications of 
the initial vacuum state or horizon-scale GR corrections 
[164] should also be tested with future measurement of 
the galaxy power spectrum. 

Upcoming observations of high redshift clusters will 
provide increased leverage on measurement of primordial 
non-Gaussianity with abundances and possibly put limits 
on any nonlinear parameter ^-t the level of a few 
tens [140]. Combining the information provided by the 
evolution of the mass function and power spectrum of 
galaxy clusters should yield constraints with a precision 
A/^L ^ 10 for a wide field survey covering half of the sky 
[215]. Alternatively, using the full covariance of cluster 
counts (which is sensitive to the non-Gaussian halo bias) 
can furnish constraints of A/^^ ~ 1— 5 for a Dark Energy 
Survey-type experiment [223, 224]. 

As emphasized in §IV however, the exact magnitude 
of the non-Gaussian bias is still uncertain partly due 
to the freedom at the definition of the halo mass and 
the uncertainty in the correspondence between simulated 
quantities and observables. Understanding this type of 
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